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ABSTRACT 

Context. Transit and radial velocity observations continuously discover an increasing number of exoplanets. However, when it comes 
to the composition of the observed planets the data are compatible with several interior structure models. Thus, a planetary parameter 
sensitive to the planet's density distribution could help constrain this large number of possible models even further. 
Aims. We aim to investigate to what extent an exoplanet's interior can be constrained in terms of core mass and envelope metallicity 
by taking the tidal Love number k 2 into account as an additional, possibly observable parameter. 

Methods. Because it is the only planet with an observationally determined k 2 , we constructed interior models for the Hot Jupiter 
exoplanet HAT-P-13b by solving the equations of hydrostatic equilibrium and mass conservation for different boundary conditions. In 
particular, we varied the surface temperature and the outer temperature profile, as well as the envelope metallicity within the widest 
possible parameter range. We also considered atmospheric conditions that are consistent with nongray atmosphere models. For all 
these models we calculated the Love number k 2 and compared it to the allowed range of k 2 values that could be obtained from 
eccentricity measurements of HAT-P-13b. 

Results. We use the example of HAT-P-13b to show the general relationships between the quantities temperature, envelope metallicity, 
core mass, and Love number of a planet. For any given k 2 value a maximum possible core mass can be determined. For HAT-P-13b 
we find Moore < 27 M e , based on the latest eccentricity measurement. We favor models that are consistent with our model atmosphere, 
which gives us the temperature of the isothermal region as ~ 2100 K. With this external boundary condition and our new k 2 -interval 
we are able to constrain both the envelope and bulk metallicity of HAT-P-13b to 1 - 1 1 times stellar metallicity and the extension of 
the isothermal layer in the planet's atmosphere to 3 - 44 bar. Assuming equilibrium tidal theory, we find lower limits on the tidal Q 
consistent with 10 3 - 10 5 . 

Conclusions. Our analysis shows that the tidal Love number k 2 is a very useful parameter for studying the interior of exoplanets. It 
allows one to place limits on the core mass and estimate the metallicity of a planet's envelope. 

Key words, planets and satellites: interiors - planets and satellites: individual: HAT-P-13b - methods: numerical 



1. Introduction 

Today more than 180 transiting exoplanets have been discov- 
ered. Within the exoplanet family planets that are found to transit 
their host star are especially important. The knowledge of mass 
and radius enables us to infer the density and bulk composition 
of a planet. Further, transit and secondary eclipse observations 
reveal information about an exoplanet's atmosphere. 

Despite all the information given, we are still far away from 
knowing the composition of the deep interior of exoplanets. 
Characteristics like the existence and mass of a potential core or 
the amount and distribution of heavy elements are quite ambigu- 
ous. Yet, this information is highly desired in order to understand 
planet formation. 

For the solar system planets, the ambiguity of interior models 
can be reduced by taking into account information from the grav- 
ity field, which is quantified by the gravitational moments J2, J a, 
and J(,. In our solar system these parameters are measureable 
and improved methods continuously increase their accuracy. 
However, for an extrasolar planet gravitational moments cannot 
be determined. Hence, we need a similar parameter that is acces- 
sible and will also provide us with information about the inte- 
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rior density distribution of the planet. This can be accomplished 
with the tidal Love numbe r k-> (lLovell91 UlG avrilov eta fl 1975b 
iGavrilov & Zharkovl[l977t IZharkov & Trubitsvnlll978h. As i t is 
equivalent to J 2 for the solar system planets dHubbardll 19841) . it 
is promising that ki will help us to further constrain the interior 
models of exoplanets. Recent studies have shown that ki is itself 
in gen eral a degenerate quantity if considered in a three-layer 
model (lKramm etal.ll2011h . 

I n this study, we fo cus on the transiting Hot Jupiter HAT-P- 
13b (Bako s~et al.l2009l) . This planet is of great interest because it 
is the only planet so far which ha s an observationally determined 
value for the Love number ko (Batyg hTet al.l 12009b . Hence, it 
gives us the chance to investigate how much information about 
the interior structure can be inferred from an actually measured 
k2 and serves as a case study in this work. We especially evaluate 
th e impact of the new eccentricity measurement for HAT-P-13b 
bv lWinn etaf]d20Toh as it gives new information about allowed 
k2 values. 

In §|2] we describe our modeling procedure, the equation of 
state (EOS) used, and our calculation of the Love number. The 
known observational constraints of the examined planet HAT-P- 
13b are given in |3] where we a l so discuss the effect of different 
eccentricity measurements, and describe the setup of our interior 
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models. The results are explained in and in $5] we compare 
with previous results. The applicability of the underlying theory 
is discussed in the Appendix. A summary is given in ^6] 



2. Methods 

2.1. Modeling procedure 

We construct planetary interior models by integrating the equa- 
tion of hydrostatic equilibrium 
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inward, where P is the pressure, r the radial coordinate, p(r) the 
mass density at r, m{r) the mass inside the radius r, and M p and 
R p the total mass and radius of the planet, respectively. We define 
the surface of the planet P p at the 1 -bar pressure level Pi = 1 bar. 
For Hot Jupiters the measure d transit radius m ay be at lower 
pressures of about ~ lOmbar dFortnev et alJl2003r) . For HAT-P- 
13b, however, the very outermost layer from 1 bar to lOmbar 
does not contribute significantly to the radius (within the error 
bars for the transit radius measurement). Hence, we start our cal- 
culation of the interior at the 1- bar pressure le vel as previously 
done for solar system planets (Guillotl fl999h . The equations 
above require the following outer boundary conditions: the ra- 
dius inside of Pi equals the planet radius, so r(Pj) 
mass inside of R p is equal to the planet's mass, so m{R p ) 
Furthermore, a planet has to obey the inner boundary condition 
m(r — 0) = 0. To fulfill this condition, we assume that a core 
exists and choose the core mass (M core > 0) such that total mass 
conservation is ensured. Equation ([TJ needs a P-p-relation (equa- 
tion of state). For details about the EOS used, see § 12.21 To obtain 
pi = p(Pi, T\) from the EOS we need the surface temperature 
T\ which we consider as a variable parameter. 

2.2. EOS 

For the interior models in this work we assume a composi- 
tion of hydrogen, helium, metals and rocks. The rocks are con- 
fined to a core and we ap ply the pressure-density relation by 
Hubb ard & Marie vl d 19891) which approximates a mixture of 
38% Si0 2 , 25% MgO, 25% FeS, and 12% FeO. For the enve- 
lope we suppose a mixture of H, He, a nd metals. For H and He 
we use the interpolated SCvH-i EOS (ISaumon et all 1 19951) . To 
represent elements heavier than He (metals) we scale the density 
of the He-EOS by a factor of four. 



2.3. Love number calculation 

For the interior models generated in this work we calculate the 
Love number Ic2 for a hydrostatic planet. This planetary property 
quantifies the deformation of the quadrupolic gravity field of the 
planet in response to an external massive body. In the cases of ex- 
oplanetary systems of int erest, the parent star of ma ss M* causes 
a tide-raising potential (Zhar kov & Trubitsvnll 19781) 
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where a is the radius of the planet's orbit, s the radial coordinate 
of a point inside the planet, 6' the angle between a planetary mass 
element at s and M, at a, and P„ are Legendre polynomials. The 
response of the planet's potential at the surface is given by 
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where R p is the radius of the planet and k„ are its tidal Love num- 
bers. Like the gravitational moments J2„ they are determined by 
the planet's internal density distribution. For a more detailed de- 
scription of the method used to calculate Love numbers and an 
analysis about some general depen dences of ki on the density 
distribution see lKramm et al.l (1201 II) and references therein. 

For the solar system giant planets the lowest order gravi- 
tational moments Ji and J\ have proven to be e xtremely valu- 
able constraints for planet modeling (se e e.g. iGuillotl fl999; 
ISaumon & Gui llot 2004; Nettelm annl201 ll) . For extrasolar plan- 
ets on the other hand, constraints derived from transit and radial 
velocity observations are often limited to the mass, radius and 
effective temperature of the planet. In a few cases information 
about the atmosphere can be obtained from spectroscopy (see e.g 
ICharbonneau et alj|2002t iPont et alJfeOQSl), which at the curren t 
level helps to validate atmosphere models dFortnev et al.l l2010). 

It can be shown that to first order in the expansion of the 
plan et's potential, th e Love number k2 is proportional to J2 (see 
e.g. lHubbardl fl984). indicating that a measurement of k2 pro- 
vides us with equivalent information like J2 for the solar system 
giants. In this way, knowledge about the Love number k2 can 
give us some insight into the planet's interior by constituting an 
additional constraint if observed. In fact, an estimation of the 
tidal Love number k2 based on observations was made for the 
Hot Jupiter HAT-P-13b, as it is part of a system in a tidal fixed 
point (see 33.11 1. The consequences of this ^-determination for 
the interior of HAT-P-13b are discussed in $4] Another possibil- 
ity to determine k2 observationally is to detect the tidally induced 
apsidal precession that i s expected to be seen in the transit light 
curves of Hot Jupiters (Ragoz zine & Wo lf 2009). 



3. Input data and model assumptions 



3.1. Observational constraints 



Discovered in 2009 bv iBakos et all (I2009). the HAT-P-13 sys- 
tem is the first exoplanetary system that contains a transiting 
planet that is accompanied by a well-characterized longer-period 
companion. The planets orbit a metal rich ([Fe/H] = +0.41 ± 
0.08 dex) G4 star with L, = 2.22 + 0.31 L . The inner planet 
HAT-P- 1 3b is a transiting Hot Jupiter which passes in front of its 
star every Pfo = 2. 916260+0. OOOOlOdays, correspondig to an or- 
bital distance of a b = 0.0427+° °^ AU. The orbit is nearly circu- 
lar with e b = 0.02 1 + O.OOSfl The radius and mass of HAT-P- 1 3b 
were determined with transit and radial velocity follow-up obser- 
vations to be M v = 0.853+°;°^ Mj, and P p = 1.281 + 0.079 Rj, 
respectively. The companion HAT-P- 13c is on a more distant 
orbit at a c = 1.188+°°^AU with a period of P c = 428.5 + 3.0 
days. As it has not been observed in transit (ISzabo et al.ll2010h 
only the minimum mass Msin/ = 15.2 ± 1.0 Mj is known. An 
important feature of the system is that the c planeQis on a highly 



(3) 



1 The value of the eccentricity is actually a matter of debate and in- 
fluences the results significantly (see j )3, 2\ . 

2 With its high mass the companion may instead be classified as a 
brown dwarf. H owever, as the defi nition of planets and brown dwarfs is 
rather unclear (Baraffe et al. 2008), we call it a planet here. 
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eccentric orbit with e c = 0.691 ± 0.018. Fo r more deta ils about 
the observed parameters of the system see lBakos et alJ (120091) . 

The high eccentricity of planet c (compared to the in- 
ner short-period planet) is one prerequisite for the theory of 
lMardUnd(l200l on the tidal evolution of a two-planet system. A 
further requirement is that the system must be co-planar. Then, 
according to Mardling's theory, the system evolves in the follow- 
ing phases: (1) The angle rj between the apsidal lines of the two 
planets circulates and the eccentricities slowly oscillate at con- 
stant amplitude accompanied by a decrease of the mean value 
of the inner planet's eccentricity. This occurs on the circular- 
ization time scale r c i rc . For HAT-P-13b r c ; rc ~ 0.05 Gyr assum- 
ing Qp ~ 10 5 and k 2 ~ 0.3. (2) rj librates around a fixed value 
r\ — 0, n. The eccentricities slowly oscillate with declining am- 
plitude but maintain the mean value of the inner planet's eccen- 
tricity. This occurs on twice the circularization time scale. At the 
end of this phase the system resides in the tidal fixed point char- 
acterized by aligned apsidal lines which then precess with the 
same rate. (3) The subsequent long term evolution proceeds on a 
time scale several orders of magnitude longer than the circular- 
ization time scale. So the tidal fixed point is actually a quasi fixed 
point. Long term tidal effects will cause a slow non-oscillatory 
decline of both eccentricities to zero unless the companion is not 
of very low mass and long period. 

iBakos et all d2009i) observed that the pericenter of HAT-P- 
13b is aligned with that of the outer planet HAT-P-13c to within 
4° + 40°. If the apsidal lines are indeed alig ned and if the system 
is co-planar, it is in the tidal fixed point. Baty gin et alJ (2009) 
(hereafter referred to as BBL) showed that under this assump- 
tion one can estimate the tidal Love number k 2 for the inner 
planet which in turn will give information about the interior of 
planet b. They first generated models that satisfy a given choice 
of the planetary mass, the radius, and an inferred planetary ef- 
fective temperature = 1649 K. For these models BBL then 
calculated k 2 from the density distribution p(r). Furthermore, 
they determined the fixed point eccentricity et, based on the con- 
straint that the orbits of both planets have an identical precession 
rate when the system is in the tidal fixed point configuration: 
Ubtouieb, h) - w CS ecuiar , where w;, tot ai is the total precession of 
the inner planet incorporating the four most significant contri- 
butions arising from secular evolution induced by planet-planet 
interaction, the tidal and rotational bulges of the planet, and gen- 
eral relativity, where ZJcsecuiar is the precessio n of the orbit o f 
planet c due to the secular evolution (see also Mardling 2007). 
BBL found that the different ^-values of the inner planet result 
in significantly different eccentricities with larger et, imply- 
ing a smaller k 2 and vice versa. They approximated this behavior 
with a fourth-order polynomial eb(k 2 ). In summary, if a planetary 
system fu l fills t he requirements for the tidal evolution theory by 
iMardlingi d2007l) and has evolved into the tidal fixed point, one 
can estimate the Love number k 2 of the inner planet by measur- 
ing the orbital parameters of the planets. It is important to note 
that the estimates for k 2 are only valid under these specific condi- 
tions, which may not apply for HAT-P-13 (see AppendixlAlfor a 
discussion of the observational evidence of the requirements for 
the tidal fixed point theory). Nevertheless, HAT-P-13b serves as 
a good example for us here to investigate what conclusions about 
the interior of an extrasolar giant planet can be drawn from an 
inferred k 2 value. 

In the case of HAT-P-13b, BBL derived from the measured 
e h = 0.021 ± 0.009 an allowed ^-interval of 0.1 16 < k 2 < 0.425 
(hereafter denoted as BBL-interval). They further concluded that 
the core mass of HAT-P-13b is 0M ffi < M c < 120M ffi . 
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Fig. 1. Relation between the fixed point eccentricity eb and the 
Love number kg. The blue l ine sh ows the fourth order polyno- 

Diff erent eccent r icity me a- 
and IWinn et al l (l2010l) are 



mial fit from IBatvgin et al.1 d2 009). 1 
surements from IBakos et alJ 0009) 



plotted in black and red, respectively, with solid lines showing 
the mean value of the measured eccentricity and dashed lines 
showing the error ba rs. The black dotted lines indicate the k 2 
interval inferred from iBatv gin et all (120091) b ased on the eccen- 
tricity measurement from Bak os et al .1 (I2009h . in the figure also 
referred to as old k 2 interval. The gray shaded area marks a re- 
gion of k 2 values not possible for models of HAT-P-13b. The 
vertical black solid line (M core = 0) shows the maximum possi- 
ble k 2 value we found based on our interior modeling (see §|4]for 
details). Com b ining the information of the e^-A^-relation from 
I Batvgin et a J (120091) . the new eccentricity measurement from 



Win n et al.l (120101) . and our interior modeling, we find a new k 2 
interval ranging from 0.265 - 0.379 (rosy shaded). 



3.2. The importance of the eccentricity measurement 

The inner planet's eccentricity is a crucial observable in our anal- 
ysis because this quantity determines the allowed k 2 interval. In 
fact, the value of th e eccentricity is somewhat uncertain. While 
IBakos et alJ d2009l) report e b = 0.021 + 0.009, a l ater s tudy of 
new radial velocity measurements by IWinn et all d2010h found 
an eccentricity of only half that value (e/, = 0.0133 ± 0.0041) . 
We favor the eccentricity measurem ent from Winn et"a"D (12010). 
because they include the data from IBakos et alJ (I2009T) and, in 
addition, use 75 new RV data points which extend the timespan 
of the data set by about 1 yr, and therefore allow a refinement 
of the orbital parameters. The consequences for the inferred k 2 
interval are illustrated in Fig.[T] 

Ba sed on the eccentricity measurement from IBakos et"ai1 
(2009), BBL found an allowed interval for k 2 of 0.1 16 - 0.425. 
This result also includes information from their interior model- 
ing, which is why the limits of the interval are not equal to the in- 
tersections of their fourth order p olynomial fit and t he error bars 
of the eccentricity measurement. lWinn et al.l (1201 Oh found a sig- 
nificantly smaller eccentricity. As the eccentricity gets smaller, 
the resulting k 2 must become larger. As a consequence, the lower 
limit of the k 2 interval is raised to 0.265. To find the upper limit, 
we make use of our interior models. We find the most homoge- 
neous model of HAT-P-13b at a maximum k 2 of 0.379. At this 
point the core mass vanishes and there are no solutions beyond 
that k 2 value (for details see ©. This upper limit is also insen- 
sitive to variations in mass and radius within the observational 
error bars (see § 14.21 ). 
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By combining the information of the fourth order poly- 
nomial fit from B BL, the new eccentricity measurement from 
IWinn et alj d2010t) . and our interior models, we derived a new 
k2 interval of 0.265 < fc 2 < 0.379. It can be expected that this 
smaller region of allowed ^-values will significantly narrow the 
allowed range of solutions for the interior of HAT-P-13b. 

3.3. Model assumptions 

For our models we use the mean values for the mass and radius 
of HAT-P- 1 3b (iBakos et alJl2009h : 

M p = 0.853Mj = 271M e , 
R p = 1.281 Rj = 14.35 R ffi . 

The uncertainties introduced by the error bars in mass and radius 
are discussed in § 14.21 We assume a two-layer structure, consist- 
ing of a rocky core and one homogeneous envelope of hydrogen, 
helium and metals. It is not necessary to consider three-layer 
models in this work because they are within the range of solu- 
tions confined by the two-layer models (see ^4]for a discussion). 
We set the helium abundance in the envelope to the solar value 
Y = 0.27 dBahcall et al.l[T995l) . The amount of heavy elements 
Z was varied from Z = 0, which gives a model with the biggest 
possible core when the other parameters are fixed, to a maximum 
value where the core vanishes. 

Due to the proximity to its star, the temperature profile of 
HAT-P- 13b is likely to be strongly influenced by the radiation of 
the star. The various processes in an exoplanet's atmosphere in- 
fluencing its temperature are very complex and beyond the scope 
of this work. As we aim to study the effects of input param- 
eters important for planet modeling on the structure and Love 
number of a planet, we vary the envelope temperature within 
the widest possible range, even though extremely cold or hot 
models may be unrealistic. For HAT-P- 13b we change the 1 kbar 
temperatures from 3800 to 9000 K. ri kbar = 3800 K gives the 
most homogeneous planet with zero core mass and models with 
rikbar > 9000 K would have fc 2 values smaller than the BBL- 
interval (see also ©. In the outermost layer of the planet from 
1 bar to 1 kbar we assume the temperature profile to be either 
adiabatic or isothermal. This is supposed to account for the un- 
certainty of how the star's radiation influences the outer envelope 
of the planet. As the planet is very close to its star, it is likely that 
there will be an isothermal layer to some extent dFortnev et alJ 
120071) . In every case, for pressures P > 1 kbar the planet is adia- 
batic. 

In addition, we compute a model series that is consisten t with 
the nongray model atmospheres from Fort nev et al.l J2007) in or- 
der to further pin down the outer boundary condition. Based on 
their pressure-temperatur e profiles of Jupiter-l ike planets around 
Sun-like stars (Fig. 3 in iFortnev et al.l 120071) . we interpolated 
a model atmosphere for HAT-P- 13b, assuming that the incom- 
ing energy flux remains constant: L«/(47ra 2 ) = const., where L* 
is the stellar luminosity and a is the planet's semi-major axis. 
Our interpolation yields T\ b ar = 2080 K as new outer boundary 
condition. It further showed that an isothermal layer may reach 
down to P a d ~ 100 bar. Our interpolated model atmosphere is 
in agreement with an atmosphere mod el specifically calcu lated 
for HAT-P- 13b, using the methods of IFortnev etal] d2007l) . The 
thickness of the isothermal layer is a variable parameter (influ- 
enced by the age of the planet). Hence, in our model series based 
on the model atmosphere, we vary P a( j as a free parameter while 
keeping T\ b al fixed at 2080 K. 




Love Number k 2 



Fig. 2. Love numbers ki and core masses of two-layer models 
of HAT-P- 13b for different 1 kbar temperatures (color coded). 
For pressures < 1 kbar the temperature profile is either adia- 
batic (colored solid) or isothermal (colored dashed). The metal- 
licity increases along a line, reaching its maximum value for 
M cole /M p = 0. The zero metallicity line for fully adiabatic mod- 
els (brown solid) separates the region of all possible planetary 
models (white) from the prohibited area (light gray). Vertical 
black dashed lines indicate the allowed interval for fc 2 -values de- 
rived by BBL. Based on that interval we find a maximum possi- 
ble core mass of 0.32 M p » 87 M®. Our new fc 2 interval is shown 
with vertical black dot-dashed lines. It lowers the uncertainty of 
the planet's core mass down to 0.1 M p « 27 M e . The other black 
lines show our models that are based on our model atmosphere 
with a fully adiabatic envelope (solid black) and P a d of 5 bar 
(dotted), 10 bar (dashed black), 50 bar (dot-dashed black), and 
72 bar (black triangle). 

Our modeling procedure c onstitutes an i mprov ement com- 
pared to the modeling done by Batvgh Vet al.l d2009l) because we 
do not approximate the core material by a constant density and 
we include the effects of different temperatures and temperature 
profiles in the atmosphere. This allows us to draw some general 
conclusions about the capability of fc 2 to constrain interior mod- 
els of extrasolar giant planets. 

4. Results 

4. 1 . Love number, metallicity, and core mass 

We present in Figs.|2]and[3]results for the calculated core mass, 
£2, and envelope metallicity. The area between an adiabatic and 
isothermal line of the same 1 kbar temperature accounts for the 
uncertainty about the outer temperature profile as the thickness 
of a potential isothermal layer is not known. Also shown are the 
models based on our model atmosphere with an isothermal layer 
of 2080 K, reaching down to 1 (fully adiabatic), 5, 10, 50, or 
72 bar, where the core disappears. In Fig. [2] the line connecting 
zero metallicity fully adiabatic models separates the region of all 
possible planetary models from the prohibited area which is not 
accessible for models of HAT-P- 13b. 

A general trend that can be seen in Fig. [2] and [3] is that the 
Love number fc 2 increases as more metals are put in the envelope. 
This is the result of fc 2 being a measure of the level of central con- 
densation of an object. Since the total mass of the planet must 
always be M p , the enrichment of the envelope with metals leads 
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Fig. 3. Love numbers ki as a function of the envelope metallicity. 
Shown are the same models as in Fig. [2] The numbers close to 
some models give the core mass in M e . The shaded areas demon- 
strate how models can be ruled out by taking into account the k^ 
intervals and assuming that the envelope metallicity is at least the 
stellar metallicity. Considering the information from the model 
atmosphere the favored models reduce to the green shaded area. 



to a decrease in core mass. A smaller core and higher metal con- 
tent in the envelope means that the planet is more homogeneous, 
which is why the Love number grows. The maximum metallicity 
is reached when the core of the planet vanishes. This marks the 
maximum possible Love number of a planet for a given temper- 
ature. 

The temperature itself also has a significant influence on the 
metallicity, core mass, and consequently Love number of the 
planet. Higher temperatures reduce the Love number h%. The 
higher the envelope temperature, the lower its density. Low den- 
sities in the envelope require a more massive core in order to 
ensure mass conservation. High envelope temperatures thus lead 
to a strong central condensation, reflecting in a small Love num- 
ber. For HAT-P-13b we find a minimum 1 kbar temperature of 
3800 K for fully adiabatic models. At this low temperature the 
zero metallicity envelope is so dense that the core is almost 
vanished (M core = 0.2 M e ). Hence, there is no enrichment of 
heavy elements possible in the envelope. At this temperature 
maximum homogeneity is reached. That translates into a max- 
imum possible Love number for HAT-P-13b of ki = 0.379, well 
below the upper limit given by the BBL-interval. At high tem- 
peratures, on the other hand, the density in the envelope is low 
enough to enable the existence of massive cores. That creates 
the possibility to enrich the envelope with metals. For instance, 
at T\ ^ai = 9000 K the core disappears at an envelope metallic- 
ity as high as 70%. The big cores of the hot models lead to very 
small Love numbers and tend to lie outside of the BBL-interval. 

We note that models with an isothermal atmosphere gener- 
ally yield smaller Love numbers than their adiabatic counter- 
parts, because they are overall hotter. 

As is obvious from Fig. [2] and [3] not all calculated models 
are placed in the BBL-interval or even in our narrower new £2 
interval, which means some of the models can be ruled out. The 
hot and metal poor models cannot fulfill the BBL-interval. Our 
new ko interval also rules out the metal rich hot models and only 
leaves us with the colder models. 

Further constraints can be made by taking into account 
model atmospheres in order to get more information about the 



planet's outer temperature profile. Our model atmosphere yields 
a 1 bar temperature of 2080 K. An isothermal layer of P a( j > 
3 bar is required to produce solutions in our new k^ interval. 

Another point we can take into account is that planets form 
out of the same cloud of gas as their host star. Hence, it is rea- 
sonable to assume that the planet's bulk metallicity is at least 
the stellar metallicity. Of course, a planet may also be enriched 
in metals due to capture of planetesimals during formation. 
Enhanced bulk metallicities for extrasolar g iant planets over 
their parent star values were indeed found (iMiller & Fortnevl 
201 1). We use this argument to also rule out interior models that 
have less than stellar metallicity (see Fig. |3). Here we can ap- 
proximate the bulk metallicity with the envelope metallicity as 
the disputable models have such small cores that most of the 
planet's metals are in the envelope. This leaves us only with 
very few models for the interior of HAT-P-13b, where we fa- 
vor the ones compatible with the model atmosphere. Combining 
all the pieces of information, we can conclude from Fig. [3] that 
HAT-P- 1 3b has an isothermal outer layer with a temperature of 
2080 K that may reach down to 3 - 44 bar. Models with isother- 
mal layers that extend to higher pressures have such small cores 
that their envelope metallicity cannot exceed stellar metallicity. 
Consequently, our models with P a[ j > 50 bar are ruled out. The 
envelope metallicity of the allowed models is 1-11 times the 
stellar metallicity. Also, the bulk metallicity cannot exceed this 
limit as M core < 27 M e . 

In order to give a better overview of the different boundary 
conditions used we plot in Fig. [4] P-T profiles of acceptable 
and some unacceptable models. Those with T\vo m < 4000 K 
and T\ kbar > 6500 K are too cold or too hot, respectively, and 
hence not shown in the plot. We find an area that roughly con- 
strains the P-T-conditions near the atmosphere that produce ac- 
ceptable models. In addition, a model also needs to have a high 
enough envelope metallicity in order to be acceptable. For exam- 
ple, our structure models compatible with our model atmosphere 
with Tibar = 2080 K and an isothermal layer down to 10 bar lie 
in the favored area marked by Fig. [4] However, models of this 
kind with Z = are not acceptable whereas models with a high 
enough envelope metallicity are, compare also to Fig. [3] 

A very important parameter describing the interior of a 
planet is its core mass. This quantity is of interest because it 
may give a hint on the formation history of the planet. Figure [2] 
shows that the relationship between a £2 value and a core mass 
is non-unique. If, for instance, &2 = 0.2 it could either be a very 
hot Tikbai = 9000 K, metal rich planet with a very small core 
or a lower temperature T\\±, m - 6500 K one with a core con- 
stituting 10% of the planet's total mass. The only information 
that can be concluded is a maximum possible core mass, given 
by the line consisting of the adiabatic zero metallicity models. 
The maximum possible core mass of HAT-P- 13b is given by 
the intersection of this line with the lower £2 boundary, yield- 
ing Moore < 0.32 M p ^ 87 M e for the BBL-interval. Taking into 
account our new &2 interval lowers the uncertainty of the planet's 
core mass down to M core < 0.1 M p « 27 M e . 

In principle, the area under each curve can be filled up with 
three-layer models (not shown in Fig. |2j. Every redistribution of 
metals in the envelope would either result in a stronger central 
condensation (smaller ki) or, if central condensation is kept con- 
stant, require a transport of material from the core to the enve- 
lope reducing the core mass. As a result of this deg eneracy intro- 
duced by a discontinuity in the envelope (see also Kramm et al. 
1201 ll) it is only possible to determine a maximum possible core 
mass. Further narrowing of the maximum core mass could also 
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Fig. 4. P-T profiles near the atmosphere of acceptable and unac- 
ceptable models. Colors and line styles are the same as in Figs. [2] 
and [3] For lines with same color and style, models with Z = 
and Z = Z max are shown. Thick lines indicate acceptable mod- 
els; compare these also to Fig. [3] The green shaded area marks 
the P-T regime for acceptable models. However, models in this 
regime are only acceptable when they also have a large enough 
amount of metals in their envelopes. 
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Fig. 5. Core mass (solid lines) and total mass of metals (dashed 
lines) as a function of envelope metallicity for different mean 
densities according to the measurement uncertainties in mass 
and radius. The numbers are the Love numbers Ic2 for zero and 
maximum possible metallicity in the envelope. Shown are only 
models compatible with our model atmosphere {T[ t, ar = 2080 K) 
for f a d = 1 bar and P^d = 50 bar. 



be achieved with a more detailed measurement of the eccentric- 
ity of HAT-P-13b as that would lead to a smaller k^ interval. 

Finally, from all the interior models we presented here, we 
favor the ones with T\ b ar = 2080 K, an isothermal outer layer of 
3 - 44 bar, and an envelope metallicity of 1 < Z/Z* < 1 1 because 
they (i) are placed in our ne w ki interval, ( i i) are consistent with 
the model atmospheres from lFortnev et al.l {2007), and (iii) have 
an envelope and bulk metallicity that is above the stellar metal- 
licity. Assuming these parameters are "correct" for HAT-P-13b, 
the planet would have a c ore mass of 27 M ffi or less. For compar- 
ison, Jupiter models from lFortnev & N ettelmannl (1201 Oh predict 



M 7 



= 15-38M ffi ,M core = 0-10M s 



and Z p = 4 - 12 Z e 



4.2. Mass and radius uncertainties 

We now investigate the effect of the measurement uncertainties 
for the mass and radius of HAT-P-13b. For that purpose, we var- 
ied the mass and radius within the Icr error-bars to obtain dif- 
ferent mean densities. As an example, we only look at models 
consistent with our model atmosphere and P a( j = 1 and 50 bar 
because all other valid models are bounded by these and react 
the same way on mass-radius variations. 

Figure [5] shows the effect on core mass and metal enrich- 
ment. The uncertainty in mean density results in an uncertainty 
of about +20 M e in core mass and mass of metals. For a given 
core mass, the envelope metallicity can vary up to +3 Z». 

There is little effect on the Love number k^. In Fig. [2] the line 
separating the prohibited area from the possible models consists 
of zero-metallicity models. Since for Z env = the Love number 
£2 only varies by about +0.003, only a minor shift of this line can 
be expected so that all the conclusions drawn for the fiducial M p 
and R p values remain valid also for other mean densities within 
the error bars. 



Table 1. Lower limits for the tidal Q value of HAT-P-13b for 
different eccentricity measurements and intrinsic temperatures 
(resulting from the allowed extension of the isothermal layer in 
the planet's atmosphere). 



Ttt [K] 


e b = 0.021 


e b = 0.0133 




(Bakos et al. 2009) 


(Winnetal. 2010) 


275 (44 bar) 


176485 


70790 


750 (3 bar) 


3 190 


1280 



4.3. Estimating the tidal Q value 

With R p = 1.281 Rj HAT-P-13b has a large radius given its mass 
and age. This implies it possesses an interior energy source, like 
many other Hot Jupiter planets. Given its eccentric orbit, this en- 
ergy source could have a strong contribution due to tidal heating. 

Based on our analysis of k2 and our interior models we found 
that the adiabatic transition is between 3 and 44 bar. We ran at- 
mosphere models to estimate what values of the intrinsic tem- 
perature (Tint) these would correspond to. At the adiabatic tran- 
sitions of 3 and 44 bar we found T- mt = 750 and 275 K, respec- 
tively. This also gives us the intrinsic luminosity L mt . Assuming 
Li D t is provided only by tidal heating , we estimated Q with the 
following relation (Mill er et al.ll2009h 



63 
T 



(GM») 



3/2 



( M,R 5 p e b 
Q 



2\ 



,-15/2 



(5) 



This gives us a lower limit on the Q value of HAT-P-13b 
since tidal heating might not be the only interior energy source. 
There could be oth er effects contributing to L mt , for instance 
Ohmic dissipation (Batvgin & Stevensonll2010t) . 

Our estimates of the lower limit on Q for the different eccen- 
tricity measurements are displayed in Tab. Q] For comparison, 
Jupiter's Q value has bee n calculated to be between 10 5 and 10 6 
dGoldreich & Soterl lT966). 
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5. Discussion 

5.1. Comparison with previous results 

In addition to an allowed ^-interval, Baty gin et al.l (|2009) also 
provided a set of possible interior models of HAT-P-13b that 
give Ic2 values within the given interval. These models range 
in core mass from to 120M ffl . However, for the given BBL- 
interval we find a maximum possible core mass that is signif- 
icantly smaller (« 87M ffi ). This differ ence is a result of dif fer- 
ent model assumptions. The models of Batygin et al.l (l2009l) all 
have cores with constant densities of 7-12 g/cm 3 (Batygin, pers. 
comm.), resembling water cores. In contrast, the models pre- 
sented inJlris_wor^_were calculated using a compressible rock 
EOS dHubbard & Marlevll 19891) . increasing the densities in the 
cores up to 12-40 g/cm 3 . As a consequence of the lower and 
constant core densities in the BBL models, a core of a specific 
mass requires a larger radius than a core of the same mass in our 
models. This results in a greater Love number than our models 
with the same core mass. That is why our models need smaller 
core masses in order to fall in the allowed ^-interval. Models of 
HAT-P-13b in this work with a core mass of 120M e would have 
a smaller £2 than the lower ^-boundary of the BBL-interval. 
H ence, the different m aximum core masses in this work and 
in Batvgin et al. (2009) arise from the different core EOS used. 
This shows the major influence of the EOS. In our example here, 
the different core EOS cause an uncertainty of about 28% in the 
maximum possible core mass. 

5.2. The eccentricity issue 

As mentioned before (see § I3.21 > the eccentricity is crucial for de- 
termining an allowed ^-interval. We used the results from BBL 
to obt ain a new hr inter val from another eccentricity measure- 
ment (IWinn et al.l 2010). Our new interval, ranging from 0.265 
to 0.379 is significantly smaller than the one previously used by 
BBL (0.116 < ki < 0.425). This allowed for a more precise 
estimate of a maximum possible core mass (M core < 27 M e ). 

It of course has to be kept in mind that the ultimate signifi- 
cance of this new eccentricity value (and deduced k^) is unclear. 
Other effects like the stellar jitter produce a large uncertainty 
on the eccentri city measurement of the inner planet, as recently 
pointed out by Payne & Fordl (1201 ll) . 

One way to help to constrain the eccentricity would be a 
detection of the occultation (secondary eclipse) of the planet 
(for instance, with Spitzer). The timing of the occultation can 
be used to place a strong constraint on ecosw, where e is 
the orbital eccentricity a nd at is the longitude of periastron 
(ICharbonneau et all 120051) . Observations to determine ecosw 
for HAT-P- 1 3b are currently being made (Harrington & Hardy, 
pers. comm.). 

6. Summary & Conclusions 

In this work, we presented new interior models of the transiting 
Hot Jupiter HAT-P- 1 3b with the aim of showing to what extent 
interior models of extrasolar giant planets can be constrained by 
using the tidal Love number k^ in addition to the known observ- 
ables mass and radius. We also varied the envelope temperature 
and metallicity in order to demonstrate the uncertainties imposed 
on the inferred interior models. 

One main result of our work is that based on the Love num- 
ber &2 one cannot draw a conclusion on the precise core mass of 
the planet. Only a maximum possible core mass can be inferred 



which is given by adiabatic zero-envelope-metallicity models. 
Taking into acc ount the allowed ^ interval (0.1 16 < ^2 < 0.425) 
derived from iBatygin et al.l (2009), we find a maximum value 
for HAT-P- Ob's core mass of 87 M m . Further radial velocity 
observations from IWinn et all (1201 Oh gave a new value for the 
planet's eccentricity, allowing us to construct a new £2 interval 
of 0.265 < £2 < 0.379. With this smaller interval we were able 
to constrain the core mass of HAT-P- 13b further down to 27 M®. 

We also showed the influence of the envelope's temperature 
and metallicity. The coldest possible model defines a maximum 
possible value for £2 as this is the most homogeneous model. For 
HAT-P- 13b we find fo.max - 0.379 . This is well below the upper 
k 2 limit from lBatvgin et al.l (120091) . 

As our analysis showed, the temperature and metallicity of 
the planet's envelope are important parameters whose knowl- 
edge would significantly increase the effectiveness of constrain- 
ing the planet's interior with the Love number A12. We there- 
fore calculated a model series based on a theoretical model at- 
mosphere from iFortnev et all (120071) . Models of this kind have 
T\ bar = 2080 K and fall in our new £2 interval if they have an 
outer isothermal layer with 3 bar < P a( j < 44 bar. With this re- 
sult we also estimated lower limits of the tidal Q value of HAT- 
P-13b. 

Based on our analysis and the applicability of the tidal fixed 
point theory we can rule out the majority of the calculated mod- 
els. Assuming our new ^-interval, we note that hot models with 
2"ikbar > 6500 K can be ruled out (compare Fig. [3}. On the other 
hand, the very cold models with T\ kbai- < 4000 K are unlikely 
because they have a metallicity that is smaller than stellar. We 
favor models with Tibar = 2080 K, an isothermal outer layer 
of 3 - 44 bar, and an envelope metallicity of 1 < Z/Z„ < 1 1 be- 
cause they (i) are placed in our ne w fc? interval, ( i i) are consistent 
with the model atmospheres from For tnev et al.l d2007l) . and (iii) 
have an envelope metallicity that is above the stellar metallic- 
ity. Assuming these conditions really apply for HAT-P- 13b the 
planet would have a core mass of 27M ffi or less. Given the cal- 
culated tidal Q, if HAT-P- 1 3b turns out to have a small core and 
an envelope enrichment similar to Jupiter, it could represent a 
Jupiter-like extrasolar planet. 

By comparing with the previous results from IBatygin et al.l 
(120091) we have seen that the core EOS affects the core mass 
and Love number, resulting into a difference of about 28% in the 
prediction of a maximum possible core mass of HAT-P- 13b. 

As pointed out in 33.21 the eccentricity of the inner planet is 
a crucial parameter because it determines the Love number k%. 
The value of the eccentricity is still unc ertain. We analyze d the 
consequences of an e b = 0.021 ± 0-009 dBakos et al.ll2009l) and 
e b = 0.0133 + 0.0041 (IWinn et alJl2010T) . 

Finally, despite the uncertainty of the inner planet's eccen- 
tricity it is still unclear whether the tidal fixed point theory from 
iMardlingl (120071) can really be applied to HAT-P- 13b. Further 
observations are necessary to prove the co-planarity and apsi- 
dal alignment configuration of the system. Also, TTV measure- 
ments could provide clues about the existence of other long- 
period companions in the system. 
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Appendix A: Applicability of the tidal fixed point 
theory 

The analysis in this wor k is based on the assumption that the 
tidal fixed point theory of Mardling (2007) can be applied in or- 
der to extract a hi value from the inner planet's eccentricity. For 
that purpose the planetary system has to fulfill the requirements 
mentioned in 93.11 Yet unproven and crucial to the evolution into 
the tidal fixed point is the co-planarity of the HAT-P-13 system. 
In fact, if the orbits of planet b and c are mutually inclined, the 
system will not evolve to the tidal fixed point. Instead it will relax 
to a limit cycle in <?/, - 77 parameter space where the eccentricity 
of the inn er planet et and the angle between the periapse lines 77 
oscillate (Mardling 2010). In that case the Love number^ of the 
inner planet can not be determined unambiguous ly. In par t icular , 
iMardlingl (1201 Ol) shows that the analysis of lBatygin et all d2009t) 
and hence the inferred ^-interval is only valid if the mutual in- 
clination of planets b and c is less than about 10°. 

A direct estimate of the mutual inclination Ai would be pos- 
sible by observing planet c in transit and measuring the Rossiter- 
McLaughlin effect. A transit of HAT-P- 13c was predicted to oc- 
cur on April 28, 2010 dWinn et al.l20 ldl). However, su ch a transit 
could not be observed. The multi-site campaign by Szabo et al. 
(2010) revealed that a transit of HAT-P-13c can be excluded 
with a significance level from 65% to 72%. Yet, a close to co- 
planar configuration cannot be ruled out. The transit probability 
of HAT-P-13c is at most 8.5% for a mutual inclination of the 
two planets of 8° and decreases to zero for inclinations < 4° 



(iBeattv & Seager! |20 ldh . Another reason why the campaign of 
Szab o et alJ (12010) could not detect a transit may be that they 
only observed for 5 days. The orbital parameters (eccentricity 
and period) of the outer planet are not well constrained due to 
the effect of the unknown stellar jitter, which translates into large 
uncertainties about the predicted transit times of planet c, mak- 
ing it necessary to observe over a wide range of nights 3 

weeks (lPavne&Fordll2011|) . 

Support for the hypothesis of co-planarity of the HAT-P-13 
system comes from the study of the Ro ssiter-McLaughlin effect 
to determine the stellar obliquity IWinn et ail (|2010) note 
that there is an indirect connection between Y*^ and Ai: a large 
Ai would lead to nodal precession of b's orbit around c's orbital 
axis, causing periodic variations in Y,^. As a consequence, it 
would be unlikely to observe a low value of Y,^ unless A; is 
also small. At any rate, only the sky-projected angle A can be 
measured (not the true obliquity whic h makes it s t ill im- 
possible to draw firm conclusions about Ai. IWinn et al I J20T0I) 
find a small A = 1.9 ± 8.6 deg, suggesting that planet b's orbital 
angular momentum vector is well-aligned with the stellar spin 
vector. This also increases the likelyhood that the orbits of plan- 
ets b and c are co-planar . 

Payn e & Fordl d20lTh also suggest a method to determine the 
mutual inclination between planets b and c with the help of tran- 
sit timing variations (TTV). In particular, they show that systems 
with Ai ~ 0° and Ai ~ 45° result in significantly different TTV 
profiles which can easily be distingui shed. In fact, H AT-P- 13b 
has been subject to a TTV campaign. iPal et alJ d201 ll) detected 
TTV in the order of 0.01 days which 'suddenly' occured during 
their observations in December 2010 and January 201 1 while in 
the first three years of observations of HAT-P- 13b there seemed 
to be a lack of TTV. If it is a periodic process, it should have 
a period of at least ~ 3 years, making it unlikely to be caused 
by HAT-P- 13c. Hence. IPdTet alJ (120 1 lh conclude that the mea- 
sured TTV may be the result o f perturbations of a nother long- 
period companion. Interestinglv. lWinn et alJ d2010l) also point to 
the possibility of an additional body d in the system as their mod- 
els of two Keplerian orbits give an unacceptable fit to the RV data 
while a model assuming a third companion with a longer orbital 
period is successful. Furthermo r e, thin king about the origin of 
the HAT-P-13 system. iMardli ng (2010) suggests that the system 
originally contained a third planet which was later scattered out 
of the system. 

Given all these issues, to our current knowledge it seems im- 
possible to say whether HAT-P- 13b is really subject t o the tidal 
evolu tion into the tidal fixed point like envisioned by Mardling 
(12001 . Further observations are necessary to pin down the mu- 
tual inclination of the planets and to prove or disprove the exis- 
tence of another companion in the system. 



3 For comparison, the solar system planets are inclined within 3.4° of 
the Earth's orbit (except for Mercury) 
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